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Motivated by the discovery of superconductivity in /3-pyrochlore oxides, we study property 
of rattling motion coupled with conduction electrons. We derive the general expression of 
the Green's function of fully anharmonic lattice vibration within the accuracy of the second 
order perturbation of electron-ion interaction by introducing self-energy, vertex-correction, and 
normalization factor for each transition. Using the expression, we discuss the characteristic 
properties of the spectral function in the entire range from weakly anharmonic potential to 
double-well case, and calculate NMR relaxation rate due to the two phonon Raman process. 

KEYWORDS: rattling, anharmonicity 
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1. Introduction 

Recently discovered /3-pyrochlore oxides, KOs 2 C>6, 
RbOs20g, and CSOS2O6, show unusual properties, which 
originate from anharmonic lattice vibration. In these 
compounds, guest alkaline ion is enclosed in an oversized 
cage made of Os-0 octahedra. As the radius of the guest 
ion gets smaller, it can move around in a larger area 
subject to weaker restoring force around its center and 
consequently shows larger atomic displacement. *' 

Resistivity of KOS2O6 shows concave-downward 
temperature dependence. 2 ' In the high temperature 
limit, electron- ion scattering dominates the resistivity, 
which is asymptotically proportional to the ion cross- 
section. By using WKB approximation, 3 ' one can show 
that thermal average of the square of the ion displace- 
ment also has concave-downward temperature depen- 
dence in the high temperature limit when the ion os- 
cillates in an anharmonic potential. NMR relaxation 
rate at the K site in KOS2O6 shows a peculiar tem- 
perature dependence, which has a low temperature peak 
and then approaches to a constant value as temperature 
is increased. 4 ' Dahm and one of the authors explained 
this phenomena by calculating NMR relaxation rate due 
to the two phonon Raman process(l/T 1 i? T) introducing 
a self-consistent quasiharmonic approximation. 5 ' They 
showed that one important aspect of the effect of anhar- 
monicity can be understood as the shift of ion spectral 
weight to higher energy with increasing temperature. 

Now, we turn our attention to the mass-enhancement 
and superconductivity as interesting phenomena con- 
cerning the coupling between electrons and anharmonic 
lattice vibration. Low-temperature electronic specific 
heat measurements show that effective mass gets pro- 
gressively large as the radius of the guest ion gets 
small. KOS2O6 shows the mass-enhancement about 
twice of that in RbOs206. 6 ' 7 ' This tendency of mass- 
enhancement depending on the guest ion radius is one 
experimental proof of the coupling of electrons with an- 



harmonic lattice vibration. On the other hand, the Ein- 
stein frequency for these compounds is not proportional 
to the inverse of the square root of the ion mass, and even 
gets smaller with decreasing ion mass. 6 ' 8 ' That is be- 
cause the restoring force gets weaker with the smaller ion 
radius in the oversized cage. Actually the band structure 
calculations confirm this tendency. 9,10 ' However, these 
compounds exhibit progressively higher T c with decreas- 
ing Einstein frequency. This unusual tendency which is 
not consistent with the well-known theory of supercon- 
ductivity is another proof of the coupling with anhar- 
monic lattice vibration, and indicates that anharmonic 
lattice vibration can possibly result in higher T c . To the 
best of the authors' knowledge, the rigorous modification 
to Eliashberg's theory 11 ' 12 ' which makes it possible to 
treat full effect of anharmonicity and reveal the relation 
between anharmonicity and electron-ion coupling does 
not exist. Then, this problem is left as an open question 
to study. As the first step toward this direction, we will 
study theoretical treatment of anharmonic ionic motion 
in this paper. 

The field theoretic method is often employed to 
study thermodynamic properties of condensed matters. 
The standard diagramatic technique is of course applica- 
ble to treat anharmonicity as the perturbation. 13 ' How- 
ever, in the extremely anharmonic case like KOS2O6 or 
tunneling state in which the harmonic part of the local 
potential is almost zero or negative, it is not reasonable to 
start from harmonic potential as the unperturbed state. 
To treat the full effect of anharmonicity, it is necessary 
for us to develop a new technique because the Wick's 
theorem is inapplicable anymore. In this paper, we use 
the set of eigenstates of effective local ion potential as the 
starting point, and treat electron- ion interaction as the 
perturbation. In contrast to the harmonic case, the tran- 
sition energy to the neighboring state is not constant, and 
transition can take place even to not neighboring states 
in an anharmonic potential. 14 ' Introducing self-energy, 
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vertex-correction, and normalization factor for each tran- 
sition with various transition energy, we will show that it 
is possible to obtain a general expression of the Green's 
function of fully anharmonic lattice vibration within the 
accuracy of the second order of the electron-ion inter- 
action. Employing that expression, we will discuss the 
characteristic properties of spectral functions of the an- 
harmonic oscillator and its application to the NMR re- 
laxation rate, and finally comment on its applicability to 
the study of mass-enhancement and superconductivity. 

2. Model and Calculation of Green's Function 

2.1 Model 

We consider the system of a local ion oscillator and 
the conducting paramagnetic electrons as described by 
the following Hamiltonian. For simplicity, we consider a 
one-dimensional oscillator. 



Hq = ^2&cl a c ka + ^2 E m \m)(m\ 



(1) 



m=0 



As usual, Cko- and cj^. are annihilation and creation 
operators of Fermions, and £t denotes a dispersion mea- 
sured from the chemical potential. E m and \m) are the 
set of eigenenergies and orthonormal cigenstates of the 
oscillator in a center-symmetric local potential V(x). 

d 2 



1 



" 2M dx 2 



+ V(x) \m) = E m \m) 



V{x) = £ 



n=l 



fan _ 
2n ' 
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(3) 



To form the bound states, &„ max must be positive. 

Let us assume that ion displacement couples linearly 
with electrons, then we write electron-ion interaction as 
follows. 



(4) 



kqcr 



where N is the total number of lattice site, and V(q) is 
determined from the inner product of the ion displace- 
ment direction and the Fourier transformation of the first 
derivative of the ion potential for the electrons. 

2. 2 Formulation of the Green's function 

We define the general retarded Green's function of 
lattice vibration as follows. 



D R (to) 



dt e iut ([x(t),x(0)}) 



(5) 



Here, x{t) is the Heisenberg representation of the ion 
displacement, [A, B] denotes the commutation relation, 
and the bra-ket (• • • } denotes the statistical average. The 
Green's function is directly related to the ion polariza- 



tion. 

On the other hand, the temperature Green's function 
is defined by employing the imaginary time Heisenberg 
operators. 



d(r) 



(T t (x(t)x(0))) 

-(x(t)x(0)) (/3>t>0) 
-(x(0)x(t)) (0>r>-/3) 



(6) 



Because of the bosonic definition and the cyclic prop- 
erty of mathematical trace, the temperature Green's 
function has the following periodic property concerning 
the imaginary time. 



d(r + f3) = d(r) 



(t<0) 



(7) 



Then, Fourier expansion is defined using the Mat- 
subara frequency u>k — 2knT, where k is an integer. 



d(iw k )= dTe^ kT d(r) 
Jo 



(8) 



By changing to the Lchmann representation, we can 
show that these two functions are related through ana- 
lytic continuation. 15 ' 

D r (lu) = d{u} + iS) {us > 0,(5 



+0) 



(9) 



This formalism above is completely the same with 
the usual harmonic phonon. 

2. 3 2nd order perturbation 

Now, we calculate d(r) for t > considering H lnt as 
the perturbation. 

T r ( e -/3(.fio + ffint) c (H + H iDt )T xc -{H a +H iat )T ^ 



d(r > 0) 



(10) 

Because the diagonal matrix elements of x are zero 
for the eigenstates in center symmetric potential, the sec- 
ond order of H mt is the lowest contribution. Introducing 
the interaction representation, 

e -(H +H iat )r = e -Hor S ^ (11) 

S(t) = 1 - f dn c H ^H int e- H ^ 
Jo 

+ f f dndr 2 c H ^H {nt c- {Ha ^- HaT ^H int c- H ^ 
Jo Jo 



(12) 



our task now is to evaluate the second order contribution. 

If the ion potential is harmonic, we can obtain arbi- 
trary order corrections employing the well-known Feyn- 
man diagrams. However, once the potential include an- 
harmonicity, we can not employ this technique anymore. 
The Feynman diagram is based on the Wick's theorem, 
which shows that the statistical average of the prod- 
uct of interaction representation operators can be de- 
composed into the product of the statistical average of 
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various combinations of two operators. To prove this 
theorem, we need the following property of the interac- 
tion representation operator, that is we can factor out 
the time dependent exponential as the classical number 



like e" 



ae 



-ujQTa a n — uqt 



a. However, we can not use 



this relation if Hq includes anharmonic terms. Consider- 
ing this fact, we now collect second order contribution 
directly without using Feynman diagram technique. 

After some algebra, we obtain a general expression 
for the second order corrections. 

fTl 

dridr 2 F(ti,t 2 ,t) 



kq 



+J2 Ei^i'/X^ 7 " 2 F (^n,r 2 ) 



kq 

JO), „ s (0) 



kq 



r r 

dridr 2 F(ti,t, r 2 ) 

T JO 



(13) 



(14) 



(15) 



kq 



m 



dridr 2 d(°)(T)d (0) ( 



X^ 0) (n -75)^(75 -Tl) 



(16) 



The free ion and electron propagators are written 
explicitly as follows for positive r, where Z i0 denotes the 
partition function of free ion part, and /(£) denotes the 
Fermi distribution function l/(e^ T + 1). 

S°\r) =-(x(t)x(0))o 

= -Z M iJ2e-^e( E -- E ^\ {m \ x \ n) f (17) 



<?f W = -(ck(r)4(Q)} = -e"^ (1 - /(&)) (18) 



9 ( uU-t) = +« + q(°) c k + q(r))o = 



/(Ck+q) (19) 



The function F(A, /x, z/) is defined as follows. 
F(A,/x,z/) = Z M X ^2 (h\x\t2) (t2\x\i3) (h\x\u) (u\x\h) 



xe 



-S i2 ) X+(E, 2 -E i3 )n+(E i3 -E H )v) 



(20) 



Now, we introduce diagrams which is different from 
the Feynman diagrams but give some hints to understand 
the physical meaning of each contribution, Eq(13)~(16). 
The rules of drawing the diagrams and correspondence 
with equations are the following. 

(1) There are two kinds of vertices, one of which is the 
external vertex with external imaginary time or 
t, and the other is the internal vertex with internal 



imaginary time t\ or r 2 . We arrange them so that 
attached imaginary time arguments proceed clock- 
wise, and connect them by wavy lines with suffixes 
of ion eigenstates. Especially for the vertex with the 
minimum imaginary time argument, we also attach 
the inverse temperature. 

(2) We connect internal vertices with two oppositely di- 
rected arrows which denote the electron-hole excita- 
tions. 

(3) We assign x-matrix elements between two ion eigen- 
states sandwiching each vertex. 

(4) Corresponding to each vertex, we multiply by the 
exponential factor in which the imaginary time is 
multiplied by the eigenenergy in smaller imaginary 
time region with negative sign, and the other one is 
multiplied with positive sign. 

(5) For two arrows, we assign free electron lines with 
certain wave vectors with imaginary time arguments 
with opposite signs. 

(6) Defining the integer I, which is the number of the 
connected diagrams, we multiply {~l) l+1 Z~^. 

(7) We multiply \V(q}\ 2 /N 2 which denotes that the in- 
teraction accompanies ±q momentum transfer at 
two internal vertices, and the factor 2 denoting the 
spin summation, and finally take summation over 
two wave vectors and ion eigenstates. 

Following this rule, we can draw diagrams as shown 
in Figs.l(a)-l(d) which correspond to Eq.(13)~(16), re- 
spectively. 




Fig. 1. Diagramatic representations of the four contributions, 
Eq.(13)~(16). The details of the notations are written in the 
main text. 
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2-4 Self-energy, vertex-correction, normalization factor 
Now we use the Fourier expansion, and arrange var- 
ious contributions by taking care of the physical mean- 
ings indicated by drawing diagrams. First of all, we con- 
sider Eq.(13) which corresponds to Fig. 1(a). If i\ = 13 
in Fig. 1(a), we can rewrite these contributions into the 
form of Fig. 2(a). The gray circle may be treated as the 
self-energy for the i\ — 14 transition. Its real part gives 
renormalization of the transition energy, and the imag- 
inary part gives its life-time. The same thing is true in 
Fig. 1(b) with Z2 = 14 as shown in Fig. 2(b). In the fol- 
lowing, we call this type of contributions "self-energy 
type". After some calculations, we can see that all the 
self-energy type contributions can be divided into two 
groups, one of which has (io>„ + E^ — E^)^ 1 as a com- 
mon factor while the other (iui n + — Ei 4 )~ 2 . 




(a) (b) 



Fig. 2. Two kinds of self-energy type contributions originating 
from Fig. 1(a) and Figl(b) respectively. 

Next we consider the following three cases, Fig. 1(a) 
with i\ 7^ Z3, Fig. 1(b) with ii ^ £4, and any type of pro- 
cesses in Fig. 1(c). We can regard these contributions as 
a modification of matrix elements of x and call "vertex- 
correction type". For example concerning Fig. 1(a), one 
can regard this contribution in two ways as shown in 
Fig. 3 (a). The external vertex with imaginary time r 
and two internal vertices combined by two wavy lines 
are united and considered as one renormalized vertex 
which is specified by the thick square in the top figure 
while it is also possible to combine the other external 
vertex with two internal vertices into one renormalized 
vertex as in the bottom figure. One can regard the top 
one as the correction to the i\-i^ transition while the 
other to i^-i^ transition. Because the product of matrix 
elements of x in the function F has i\-i^ permutation 
symmetry, it is always possible to re-interpret 13-14 type 
contributions as type one by permuting i\ and 13. 
In the course of calculation, we decompose every term 
of Fourier expansion coming from Eq.(13) into partial 
fractions and divide into two groups, one of which has 
{\ijj n + Ei x — Ei^)^ 1 as a common factor and the other 
(iuj n + Ei 3 — Ei A ) . About the latter one, we permute 
13 and i\ so that all the terms have the common factor 



(ko n + Ei 1 — _Ei 4 ) _1 which correspond to the correction 
to the transition. Similarly, Fig. 3(b) and 3(c) indi- 
cate two ways of consideration of Fig. 1(b) with 12 7^ «4 
and Fig. 1(c) respectively. In both cases, we follow the 
similar procedure mentioned above. When we decompose 
the Fourier expansion of Eq.(14) and Eq.(15) into partial 
fractions, we notice that all terms can be divided into two 
groups in each case. From Eq.(14), we obtain two groups 
with two kinds of common factor (ioj n + — E 1 ^) -1 or 
(iu) n + E il — Eiz)^ 1 . About the latter one, we permute 
i2 and ii considering the permutation symmetry of the 
product of matrix elements of x in the function F. On the 
other hand from Eq.(15), we will obtain two groups with 
(\<jj n + Ei t — £"i 4 ) _1 or (iuj n + Ei 2 — Ei 3 )~ x as a common 
factor. We permute i^ and i%, i^ and i^ simultaneously 
concerning the latter one. Eventually, we can arrange 
"vertex-correction type" so that all the terms have the 
common factor (iw n + Ei x — E 1 ^) -1 which correspond to 
the correction to ii-ii transition. 

The remaining part of Eq.(16) which is represented 
by the "disconnected diagram" Fig. 1(d) will contribute 
just as the total normalization of S ^ . 



*4 44 *4 




Fig. 3. Three kinds of vertex-correction type contributions. See 
the main text for more details. 

Now, we define the self-energy II, vertex correction 
Sx, and normalization factor A for each transition. 

ilii 

^ I <zi |a;|z 4 ) | 2 (1 + A(?i,i 4 )) + 5x{iuj n , h, ij) 
\uj n + E it - E u - n(iw n , h,u) 

Within the accuracy of the second order, this form 
is equal to the following Taylor expansion. 

d(iw n ) = d<°>(iw B ) 
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W EOT*** - e-^) W t l-^^y ) (22) 

+^o i E(^ i - e "^) J x( r ; ; ,Z1 ' f ( 23 ) 



e -/3E I4 j |<^l|^N 4 )| 2 A(ii,M) 

iw„ + E il - E iA 



(24) 



We can derive the expressions of II, Sx, and A by the 
following procedure. About the self-energy H(iuj n , i\ 1 14), 
we collect the terms with the common factor (iu> n + 
E il — E i4 )~ 2 in the self-energy type corrections and 
take them equal to Eq.(22). Next, we collect all the 
terms of vertex-correction type as the corresponding 
part of Eq.(23), then we can obtain the expression of 
5x(kj n , ii, 14). Finally about the normalization factor 
we collect disconnected diagram contribution 
and (icj n + E il — Ei^)" 1 liner terms of self-energy type 
corrections and take them equal to Eq.(24). The final 
results are shown in the Appendix explicitly. 

There are some notes worth paying attention to con- 
cerning the normalization factor. In the harmonic case, 
we can show that the disconnected diagram contribution 
and linear terms of the self-energy type corrections 
exactly cancel out when we take summation over ion 
eigenstates. This cancellation corresponds to the linked 
cluster theorem in the Fcynman diagramatic technique. 
Additionally, this cancellation also occurs in the anhar- 
monic case if we take the low-temperature limit. Con- 
sidering this fact, it is reasonable to introduce normal- 
ization factor in the numerator in Eq.(21) and include 
these terms in A simultaneously. About the nonlinear 
terms with the common factor (iu n +E il — E^)^ 1 which 
comes from the self-energy type corrections originally, it 
may be ambiguous whether they should be treated as 
the self-energy or normalization factor although these 
two ways are equivalent in the second order perturba- 
tion. Actually, those terms are ui independent and gives 
only an unimportant corrections unless the perturbation 
fails. Then this ambiguity is not a serious problem. 

At the end of this section, we would like to make 
some comments of the preceding study by Dahm and one 
of the authors. 17 > In the paper, they treated the square- 
well potential as an example of extremely anharmonic 
case, and drew spectrums introducing self-energy phe- 
nomcnologically which is common to all the transition 
energies. They succeeded to describe the spectral weight 
shift arising from the anharmonicity, but the obtained 
spectral function includes unphysical zero points which 
correspond to the singularities between various transi- 
tion energies in the free propagator. By introducing the 
self-energy, vertex-correction, and normalization factor 
for "each transition" as discussed above, we can avoid 
the unphysical phenomena. 



3. Spectral Function and NMR Relaxation Rate 

We draw spectrums for several typical cases using 
the retarded Green's function which is obtained by the 
analytic continuation. 

1 



%) = -; 



r ImDji(w) 



(25) 



7r|(0Ml)P 

Introducing appropriate scale of length x = dx and 
ion energy scale E, we write the local ion Hamiltonian 
as follows. 



• ax 2 | /„•' 



(26) 



We take the constant E = 10 as an typical energy scale 
of Einstein frequency with milli-electron volt in mind as 
the units, and change two parameters a and 6. If a > 0, 
we keep the relation a+b = 1 so that the potential energy 
ax 2 + bx 4 in the units takes a value of unity at x = 1. 
If a < 0, we change only a keeping 6=1. Especially if 
a = 1 and 6 = 0, the Einstein frequency of this harmonic 
oscillator is 20. 

Later on, we neglect the momentum dependence of 
|V(q)| and replace by a constant value g for brevity 
of calculation. We also assume the constant density of 
states. Except for nonlinear terms in the normaliza- 
tion factor, we can perform the summation over wave- 
vectors analytically at absolute zero temperature. We 
neglect temperature dependence of those terms and em- 
ploy absolute-zero analytic expression because the en- 
ergy scale of electrons is much larger than that of the 
ion. We choose parameters, half of the bandwidth (Fermi 
energy measured from the band-bottom) = 100 keeping 
large mass-enhancement in mind, and coupling constant 
and scale of length dg = 10. In the following calculation, 
we will neglect nonlinear terms in the normalization 
factor which will give only a few percent correction with 
these parameters since it is time-consuming to perform 
numerically the integration involved. 

Close to the harmonic limit like Fig. 4, almost de- 
generate peaks sit and the spectral weight shifts to 
higher energies with increasing temperature. This shift 
is qualitatively consistent with the description of the 
self-consistent quasiharmonic approximation. 5 ) In this 
limit, many peaks are united and it appears that 
the damping is enhanced effectively. As anharmonicity 
increases(Fig.5^7), excitation energy to the neighbor- 
ing state becomes large progressively as the energy of 
the initial state increases. Then, we observe the sequence 
of peaks of which the weight shifts to higher energy 
monotonously with increasing temperature. As the co- 
efficient of the square term of x become large to negative 
value(Fig.8— »11), quasi degeneracy between the bonding 
and anti-bonding states appears. In this range, we can 
divide transitions into two groups, one of which expresses 
the transition between quasi-degenerated states lying in 
lower energy range, and the other expresses transitions 



5/?? 



J. Phys. Soc. Jpn. 



Full Paper 



s. 
< 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




2.4 
2 
1.6 



18 20 22 



[0-1] 




-2 2 
Displacement(d) 



3^ 
< 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




2.4 

2 
1.6 



24 26 



1] 



T=2 
T=5 
T=10 
T=20 
T=40 




IV2] 



-2 2 
Displacement(d) 

[2-3] 



10 



20 



30 



40 



50 



10 



20 



30 



40 



50 



Fig. 4. The spectral function of a local ion oscillator(Eq.(26)) for 
the parameters (a,b)=(0.98,0.02). The inset at the top comple- 
ment the lowest energy excitation peak, and the other one shows 
the potential and eigenenergies. The eigenstate suffix runs from 
0. For each peak, two numbers are attached to specify transition. 
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Fig. 5. (a,b)=(0.9,0.1) The notations used in the figure are the 
same as Fig. 4 



between non-degenerated states lying in higher energy 
range. An interesting case is the shoulder at lj ~ 30 
seen in Fig. 10. This characteristic spectral shape is due 
to the crossover of the above mentioned two sequences. 
Between a = —3 and a = —4, the transition energy de- 
noted by [2,3] gets smaller than that of [1,2], then non- 
monotonous tendency of the peak height which corre- 
sponds to the Boltzmann weight of the initial states ap- 
pears. In the off-center limit (a — > — oo), the ion eigen- 
states approaches the linear combination of two off-center 
oscillations. Peaks will be consolidated to two peaks of 
which the lower one approaches zero energy and the other 
denotes the energy of the off-center oscillation. 

NMR relaxation rate due to the two-phonon Raman 



Fig. 6. (a,b)=(0.5,0.5) The notations used in the figure are the 
same as Fig. 4 



J3 
< 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




2.4 
2 
1.6 



[C 1] 



26 28 




-2 2 
Displacement(d) 



T=2 
T=5 
T=10 
T=20 
T=40 



[1:2] 



[2-3] 



10 



20 



30 



40 



50 



Fig. 7. (a,b)=(0,l) The notations used in the figure are the 
same as Fig. 4 



process is expressed as follows. 5 ) 
1 



t r 



cx / duj A 2 (uj) [n{us) + 1] n(u>) 



(27) 



where n(u>) denotes the Bose distribution function. The 
plot of l/Tf-T for various type of potential is shown 
in Fig. 12 and 13. In the harmonic limit like (a,b) — 
(0.98,0.02), 1/Tf-T increases monotonously. Then, with 
increasing anharmonicity, we see that the low-energy 
peak which is observed in the experiment 4 -' develops. 
That change is consistent with the result of the self- 
consistent quasi-harmonic approximation. 5 ) 

The behavior of 1 /TfT with increasing anharmonic- 
ity is explained in relation to the spectral function as 
follows. If anharmonicity is strong enough, the spectral 
weight shift to higher energy takes place faster than the 
gain of Bose function weight. Then low energy peak ap- 
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Fig. 8. (a,b)=(-l,l) The notations used in the figure are the 
same as Fig. 4 



Fig. 10. (a,b)=(-3,l) The notations used in the figure are the 
same as Fig. 4 
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Fig. 9. (a,b)= 
same as Fig. 



(-2,1) The notations used in the figure are the 



Fig. 11. (a,b)=(-4,l) The notations used in the figure are the 
same as Fig. 4 



pears as shown in Fig. 12. As the low frequency peak in 
the spectrum grows towards the off-center limit, the peak 
in the NMR relaxation rate gets sharpened as shown in 
Fig.13. 

4. Summary and Discussion 

In this paper, we have derived the general expres- 
sion of the ion Green's function in the second order per- 
turbation of electron-ion coupling. Explicit forms of the 
self-energy, vertex-correction, and normalization factor 
are obtained independent from the strength of the an- 
harmonicity. Although we treated the oscillator moving 
in one dimension, it is straightforward to generalize the 
results to a three dimensional oscillator. In this sense, 
the present theoretical scheme may be useful for a wide 
class of problems concerned with anharmonic lattice vi- 
bration. 



Using the expression, we have discussed characteris- 
tic properties of ion spectral functions in the entire range 
of anharmonicity, from almost harmonic case to the well- 
separated double- well potential. These spectrums are ob- 
servable by THz absorption and we expect that experi- 
mental results 18 ^ can be understood qualitatively as the 
combination of characteristics discussed in the previous 
section for some typical cases. Additionally, we calcu- 
lated NMR relaxation rate due to the two phonon Ra- 
man process, and reproduced the low temperature peak. 
The spectrum of the anharmonic lattice vibration is nat- 
urally reflected in other experiments. The present work 
will provide a basis of analysis of ultrasonic measure- 
ments 19-21 - 1 and Raman scattering experiments 22-24 -' on 
/3-pyrochlore, skutterudite, and clathrate compounds. 

The higher order perturbation is possible at least in 
principle using continued fractions, but it will be cum- 
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Fig. 12. NMR relaxation rate due to the two phonon Raman pro- 
cess at the ion site(Eq. (26)) with a > 0. 



0.18 
0.16 
0.14 
0.12 
0.1 



0.08 
0.06 
0.04 
0.02 




(a,b)=(0,1) 
(a,b)=(-1,1) 
(a,b)=(-2,1) 
(a,b)=(-3,l) 
(a,b)=(-4,1) 



10 



15 
T 



20 



25 



30 



Fig. 13. NMR relaxation rate due to the two phonon Raman pro- 
cess at the ion site(Eq. (26)) with a < 0. 



bersome and laborious task. About the electronic proper- 
ties, it may be straightforward to apply similar procedure 
based on the second order perturbation expansion. The 
merit of the method in this paper is that it allows us to 
consider full effect of anharmonicity as we have stressed 
already. In the weak coupling region, it may be possible 
to construct a microscopic theory which will reveal the 
relation between anharmonic lattice vibration and elec- 
tron mass-enhancement, and furthermore superconduc- 
tivity eventually. However, we will leave these interesting 
problems as future study. 
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Appendix: Explicit expressions of Self-energy, 
vertex correction, and normalization 
factor 

We obtain the expression of self-energy, vertex cor- 
rection, and normalization factor as Eq.(A-l),(A-2),(A-3) 
by following the procedure written in the Sec. 2. 4. 
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